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Stochastic Structural Stability Theory (SSST) provides an autonomous, deterministic, nonlinear 
dynamical system for evolving the statistical mean state of a turbulent system. In this work SSST 
is applied to the problem of understanding the formation of the roll/streak structures that arise 
from free-stream turbulence (FST) and are associated with bypass transition in boundary layers. 
Roll structures in the cross-stream/spanwise plane and associated streamwise streaks are shown to 
arise as a linear instability of interaction between the FST and the mean flow. In this interaction 
incoherent Reynolds stresses arising from FST are organized by perturbation streamwise streaks to 
coherently force perturbation rolls giving rise to an amplification of the streamwise streak pertur- 
bation and through this feedback to an instability of the combined roll/streak/turbulence complex. 
The dominant turbulent perturbation structures involved in supporting the roll/streak/turbulence 
complex instability are non-normal optimal perturbations with the form of oblique waves. The 
cooperative linear instability giving rise to the roll/streak structure arises at a bifurcation in the 
parameter of STM excitation parameter. This structural instability eventually equilibrates nonlin- 
early at finite amplitude and although the resulting statistical equilibrium streamwise streaks are 
inflectional the associated flows are stable. Formation and equilibration of the roll/streak struc- 
ture by this mechanism can be traced to the non-normality which underlies interaction between 
perturbations and mean flows in modally stable systems. 



I. INTRODUCTION 



The physical mechanism of turbulence in shear flows is not yet comprehensively understood despite many recent 
advances in experiment, simulation and theory. The problem of shear flow turbulence can be divided into two 
components: transition from the laminar to the turbulent state and maintenance of the turbulent state. The transition 
problem results from the lack of an inflection in the velocity profiles of most boundary layer flows. Inflections are 
associated by the Rayleigh theorem with existence of robust instabilities that continue in viscous flows from the 
inflectional instability of the same velocity profile in an inviscid flow. The problem of observed robust disturbance 
growth in perturbation stable shear flows was solved when it was recognized that the non-normality of the underlying 
linear dynamics of shear flows allows perturbation growth in the absence of exponential instability. The concept of 
transient growth in shear flow has roots in the classical work of Kelvin and Orr [1, 2] who used analytical solutions 
of perturbation dynamics in idealized shear flows to provide example solutions demonstrating transient growth. This 
early work remained obscure presumably due to lack of a convincing physical application. Transient growth concepts 
were first applied in the modern context of linear operator non-normality to understanding the three dimensional 
baroclinic turbulence in the midlatitude atmospheric jet after comprehensive data collected beginning in the middle 
of the last century for the purpose of weather forecast provided convincing evidence that turbulence in the midlatitude 
jet stream was maintained by growth processes unrelated to exponential instability [3-5]. In the context of laboratory 
shear flows although three dimensional perturbations in the form of a roll/streak structure were observed in boundary 
layers (Townsend, 1956; Kline et al., 1967; Blackwelder and Eckclmann, 1979; Robinson, 1991) and related to the 
nonmodal lift-up growth mechanism [6] similarly comprehensive observational evidence for the mechanism of nonmodal 
growth in boundary layer flows awaited the advent of direct numerical simulation (DNS) at Reynolds numbers 0(1000), 
for which turbulence is maintained in shear flow, and particle image velocimetry (PIV) of turbulent laboratory shear 
flows. The methods of non-normal operator analysis and optimal perturbation theory were first applied in the 
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context of laboratory shear flows to two dimensional disturbances [7]. It was believed at the time that secondary 
instability of finite amplitude two dimensional equilibria were the mechanism of transition [8, 9] and it was shown 
that these unstable two dimensional nonlinear finite amplitude equilibria could be readily excited by even very small 
optimal initial perturbations[10]. However, it became increasingly apparent from observation and simulation that 
the finite amplitude structures associated with transition are three dimensional and analysis of three dimensional 
optimal perturbation growth followed [11-17]. These analyses revealed that the optimally growing three dimensional 
structure is associated with cross-stream/spanwise rolls and associated strcamwise streaks and is related to the linear 
lift up mechanism. The remarkable convected coordinate solutions for perturbation growth in unbounded shear flow 
[1] allow closed form solution for the scale independent structures producing optimal growth in three dimensional 
shear flow [13, 14]. These closed form optimal solutions in unbounded shear flow confirm the result found numerically 
in bounded shear flows that for sufficiently long optimizing times streamwise rolls produce optimal energy growth 
while for short optimizing times the optimal perturbations are oblique wave structures that synergistically exploit 
both the two dimensional shear and the three dimensional lift up mechanisms producing vortex cores oriented at an 
angle of approximately 60 degrees from the spanwise direction. And indeed, the roll/streak and oblique accompanying 
structure complex that is predicted to produce optimal growth by analysis of non-normal perturbation dynamics of 
shear flows has been convincingly seen in both observations and simulations [18-21] and shown to be essentially related 
to the non-normality of shear flow dynamics [22, 23]. 

Although the mechanism of non-normal growth has been clarified, and its importance in bypass transition and 
maintenance of the turbulent state is widely if not universally accepted, the route by which non-normality leads to 
formation of the roll/streak structure and the part played by this coherent structure, its nonlinear equilibration and 
its stability, in the transition to turbulence and maintenance of the turbulent state remains to be determined. 

The most direct mechanism exploiting non- normality to form roll/streak structures is introduction of an optimal 
perturbation into the flow, perhaps by using a trip or other device [11, 15, 16]. A related approach is to stochastically 
force the flow, with stochastic forcing regarded as modeling surface roughness or free stream turbulence (FST). 
This mechanism can be analyzed using stochastic turbulence modeling (STM) [24-31]. Because of the non-normal 
nature of perturbation growth in shear flow, stochastic turbulence models are closely related to optimal perturbation 
dynamics. In conventional stochastic turbulence models the roll/streak structure is envisioned to arise from chance 
occurrence of optimal or near optimal perturbations in the stochastic forcing[29, 32-34]. These mechanisms exploit 
the linear non-normal growth process directly. However, the ubiquity of streak formation suggests, as argued by 
Shoppa and Hussein [19], that some form of instability process underlies the formation of streaks, that this instability 
involves an intrinsic association between the roll/streak structure and associated oblique waves and vortices, and that 
this three dimensional instability must differ qualitatively from the familiar laminar shear flow instability. Indeed, 
from a comparison of experiment with simulation Anderson ct al. [35] conclude that the evidence "..corresponds 
to some fundamental mode triggered in the flat-plate boundary layer when subjected to high enough levels of free- 
stream turbulence..". Previously proposed exponential instability mechanisms include centrifugal instability [36] and 
the Craik-Leibovich instability [37]. Proposed algebraic growth mechanisms involve a streamwise average torque 
produced by interaction of discrete oblique waves [38, 39]. 

All these streak growth mechanisms rely on the dominant non-normal process in shear flow which is lift up of mean 
streamwise velocity by perturbation cross-stream velocity, while differing in the manner in which this cross-stream 
velocity arises. 

The cross-stream/spanwise roll structure provides a powerful mechanism for forming streamwise streaks in shear 
flows whether episodically forced by an initial condition or continuously forced by an oblique wave structure. However, 
in the absence of feedback from the streak back to the roll this powerful streak formation mechanism does not result 
in instability although because of the large streak growth produced by a cross-stream/spanwise roll perturbation, 
placing even a very weak coupling of the streak to the roll, such as by a small spanwise frame rotation, results in 
destabilization [40, 41]. Turbulent Reynolds stress provides an alternative mechanism for producing the feedback 
between the streak and roll needed to destabilize the roll/streak structure. Indeed, if we observe a turbulent shear 
flow in the cross-stream/spanwise plane at a fixed streamwise location we see that at any instant there is a substantial 
torque from Reynolds stress divergence forcing cross-stream/spanwise rolls. The problem is that this torque is not 
systematic and so it vanishes in temporal or streamwise average. However, in the presence of a perturbation streak the 
symmetry in the spanwise direction is broken and the torque from Reynolds stress divergence can become organized 
to produce the positive feedback between the streak and roll required to destabilize this structure by continuously and 
coherently exploiting the powerful non- normal roll/streak amplification mechanism. The existence of this mechanism 
for destabilizing the roll/streak structure in turbulence makes it likely that some dynamical perturbation complex 
exists to exploit it. In this work we prove by construction that this is so by deriving a system of equations eigenanalysis 
of which reveals the unstable roll/streak/turbulence structure that is responsible for destabilizing turbulent shear flow 
to streak formation. This emergent instability can be understood as a synthesis of the streak formation mechanisms 
described above in which FST, rather than itself constituting the cross-stream velocity linearly forcing streak growth, 
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instead is organized by the perturbation streak into oblique waves that quadratically force the cross-stream/spanwise 
roll by inducing a Reynolds stress torque linearly proportional to streak amplitude thereby producing an emergent 
exponential instability of the combined roll/streak/turbulcncc complex. As a boundary layer develops in the presence 
of FST this exponential instability of organization of the roll/streak complex is the first non- viscous instability to 
occur. We remark that this coupling in the highly non-normal boundary layer shear flow of the roll and streak 
components by turbulent stresses resulting in destabilization of the roll/streak structure is analogous to the coupling 
of the torroidal and poloidal components of the magnetic field by the turbulent a effect which destabilizes the induction 
equation in the magnetic dynamo problem [42]. 

The challenge is to find a method of stability analysis analogous to the method of modes for laminar flow instability, 
but applicable to the emergent turbulence/mean flow interaction instability. Specifically, analysis of the cooperative 
instability of the roll/streak/turbulence complex requires constructing a dynamical system for evolving the consistent 
statistical mean of a turbulent state. What is required is a physically correct while at the same time a computationally 
accessible approximation to the evolution of the trajectory of the probability density function (pdf) of a turbulent 
system in phase space. We refer to this dynamical system as the stochastic structural stability theory (SSST) system. 
The approximate pdf of the turbulent state, which is carried forward in time by the SSST system dynamics, often 
evolves to a fixed point corresponding to a statistically steady turbulent state. This method for analyzing the dynamics 
of the pdf of a turbulent system was developed to study the phenomenon of spontaneous jet formation at global scale 
in planetary atmospheres [41, 43, 44] and has also been applied to the problem of spontaneous jet formation from 
drift wave turbulence in magnetic fusion devices[45]. In SSST the turbulence is simulated using a STM in which 
intrinsic excitation of perturbations by nonlinear scattering and extrinsic excitation by FST are parameterized as 
stochastic [46-49]. The STM provides an evolution equation for the quadratic statistics of the turbulent eddy field 
associated with a mean flow. In the STM the eddy field is expressed in terms of a covariance matrix from which the 
Gaussian probability density function approximation for the turbulence variance and quadratic fluxes can be obtained. 
Coupling a time dependent STM to an evolution equation for the streamwise mean roll/streak/shear complex produces 
a nonlinear dynamical system for the co-evolution of the roll/streak/shear and the self-consistent quadratic statistics 
of its associated turbulence: this is the SSST system. 

The SSST equations incorporate a stochastic turbulence model but these equations are themselves deterministic 
and autonomous with dependent variables the streamwise mean roll/streak/shear complex and the streamwise mean 
covariance of the turbulence. The perspective on shear flow stability provided by these equations differs from the more 
familiar perspective based on perturbation stability of stationary laminar flows. In fact, the primary perturbation 
instability in SSST has no counterpart in the stability theory of laminar flow; it is rather a cooperative instability 
in which the evolving roll/streak/shear complex organizes the background turbulence covariance to produce flux 
divergences configured to amplify the roll/streak/shear complex leading to an emergent coupled roll/streak/shear plus 
turbulence instability that does not involve perturbation instability of the streak. The SSST equations approximate 
the nonlinear streamwise mean dynamics of the coupled roll/streak/shear plus turbulence complex and this system in 
many cases supports equilibration of the emergent roll/streak/shear complex and its consistent turbulence field at finite 
amplitude. Equilibrium between a mean flow and its associated field of turbulence requires that the momentum flux 
divergence arising from the turbulent and mean velocities produce a stationary state of balance with the streamwise 
mean flow forcing and dissipation. The remarkable fact is that the turbulence, which depends on the roll/streak/shear 
complex, and the roll/strcak/shcar complex, which depends on the turbulence, quite generally adjust to produce such 
balanced states. 

While the solution trajectory of an initially unstable SSST state generally converges to a fixed point representing 
a state of balance among the mean flow forcing and advection, the turbulent Reynolds stress divergence, and the 
damping; these finite amplitude equilibria may lose structural stability as a function of the STM excitation parameter 
and this instability then leads either to another equilibrium or to a time dependent limit cycle or chaotic solution 
[41, 43, 44, 50, 51]. Chaotic trajectories of the SSST system correspond not to the familiar chaos of an individual 
turbulent state trajectory but rather to chaos of the ensemble mean turbulent state. An example of this kind of chaos 
is the irregular fluctuation of the mean flow/turbulence complex seen in drift wave turbulence [52]. 

In this work we concentrate on the emergence of roll/streak structures in boundary layers as an instability of 
interaction between FST and the streamwise mean flow and on the mechanism by which this instability equilibrates at 
finite amplitude to maintain a statistically stable mean roll/streak/turbulence structure. Stable roll/streak structures 
are initiated as an instability at a minimum intensity of FST and are maintained as statistically stationary structures 
over an interval of STM excitation parameter. The time dependent state initiated at higher levels of FST will be 
reported on elsewhere. 
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II. REVIEW OF STREAMWISE VORTEX FORCING BY OBLIQUE WAVES 



The emergent SSST streak instability mechanism can be understood as a synthesis of mechanisms previously inves- 
tigated in the study of streak formation in boundary layers and we begin by reviewing the mechanism of strcamwise 
vortex forcing by oblique waves which has been used [38, 39, 53] to explain the streamwise mean vortex circulations 
observed by Klcbanoff ct al. [18]. These works ascribe streak formation to interaction of a specific mixture of Tollmicn- 
Schlichting (TS) waves and oblique perturbations chosen to produce through quadratic advcctive interaction a zero 
streamwise wavenumber component. Consider a fluctuating eddy field and its associated Reynolds stresses UiUj, where 
Ui is one of the three components of the eddy velocity field, u in the streamwise, x, direction; v in the cross-stream, 
y, direction; and w in the spanwise, z, direction. Here and in the sequel an overline denotes the streamwise average. 
Divergence of the Reynolds stress induces a streamwise mean force: 

_ dpuTuJ 
dxj 

The mean spanwise force is 

dpvw dpw 2 
dy dz 

and the mean cross-stream force is: 

dpvw dpv 2 
dz dy 

Consequently the streamwise component of the mean torque G = V x F is: 



Gx ^ dF^^dFy 

dy dz 

dz 2 dy 2 ) ^ dydz^\ 

If the eddy field is spanwise homogeneous no streamwise mean torque is present. However, it is a remarkable fact 
that quite generally in three dimensional flows when spanwise symmetry is broken mean torques arise that maintain 
streamwise mean vortices. 

As an example it is instructive to review the case of an unbounded domain with a flow field consisting of two 
oblique structures in the x, z plane with wave vector (k, m) inclined at an angle = ± tan -1 (m/k) to the streamwise 
direction: 

u = u cosmze lkx+dy , v = v sinmze lkx+Uy , w = w cosmze lkx+dy , (1) 
With velocity components satisfying non-divergence: 

ikuo + Uvq + mvjQ = . 

The Reynolds stresses for this flow are: 



vw = —^R,(voWq) sin2mz , v 2 = — | | 2 cos 2 mz , w 2 = —\wq\ 2 sin 2 mz . 

The spanwise dependence of the Reynolds stresses implies a streamwise mean force F y (z) in the cross-stream plane 
inducing streamwise mean toque: 

G x = — pm 2 ^{vQWQ) sva.2mz . 

Note that this torque has double the spanwise wavenumber of the perturbation field. There is no streamwise torque 
induced by a single oblique wave when m = 0. Perhaps less obvious is that no Reynolds stress torques are induced 
by streamwise rolls with k = because continuity requires vq and wq be in quadrature and therefore ^.(vqWq) = 0. 
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The vanishing of the torque in these two limiting cases is made clear by using the continuity equation to write 

JzTIl? 

G x = —p — - — SKuqWq) sin2mz , 

it is also immediate from this expression that for given square wavenumbcr k 2 + m 2 the torque is maximized for 
oblique waves with 9 = tan _1 (m/fc) = 54.7°. 

For unbounded constant shear flow, U = ay, plane wave solutions in closed form exist [1]. For such a flow it can 
be easily shown that the Reynolds stress for an initial eddy field in the form of the oblique perturbations in the form 
(1) ^: " 

,^ 1 k 2 + l 2 + m 2 + l 2 _„ m „, . n 

VW{t) = 4 k 2 + m* + (l-akt)* e 9(t) ^o)^2m Z , 

where e~ 9 ^> is the decay due to dissipation. The time dependent Reynolds stress produced by an oblique wave 
in constant shear flow produces torque with constant spatial structure and as a result coherent mean streamwise 
torque over its lifecycle. Also it can be shown using this closed form solution [f 3] that energy growth over short time 
intervals is maximized for oblique perturbations with = 63° which is very close to the orientation of the oblique 
wave maximizing the induced torque. It is consistent then that near <d — 63° lies the oblique waves that producing 
the greatest mean streamwise torques when integrated over their life cycle. 

However, the fact that a specific pair of plane oblique waves force streamwise vortices does not explain the presence of 
streamwise vortices when the mean flow and the perturbation field statistics are spanwise uniform. Remarkably, while 
no mean streamwise torque arises from turbulence in a fluid with only a constant shear in the cross-stream direction, 
the slightest spanwise structure in the streamwise flow induces coherent streamwise mean torques. An example of this 
phenomenon obtained by imposing the mean streamwise perturbation test function SU = e cosijry / (2L y )) cos(2irz/L z ) 
on a Couette flow is shown in Fig.l. This result was obtained using the STM described in section III by first calculating 
the equilibrium perturbation velocity covariancc maintained by stochastic excitation in the perturbed streamwise mean 
flow and then obtaining from this perturbation velocity covariance the resulting Reynolds stress induced torque. This 
coherent torque arises because the eddy field is modified by the spanwise variation of the streamwise streak. This 
induced torque in turn modifies the mean flow which by interaction with the eddy field produces a modified torque. 
Most mean flow perturbations organize torques that do not exactly amplify that mean flow perturbation, as is the case 
for the perturbation in Fig. 1. However, exponential growth of streamwise mean flows and associated eddy fields would 
result, at least for sufficiently small perturbations, if the mean flow perturbation were to organize precisely the eddy 
field required for its amplification. We show constructively below that these cooperative instability structures arise 
quite generally by using an eigenanalysis to solve for them. Boundary layer flows in general support this structural 
instability in the presence of sufficient FST. For a range of FST intensities this structural instability equilibrates at 
finite amplitude producing stable but strongly inflected spanwise streaks. 



III. FORMULATION OF THE COMPOSITE ROLL/STREAK/TURBULENCE DYNAMICS 

We decompose the velocity fields into streamwise mean components (indicated uppercase) and perturbations (in- 
dicated lowercase) so that the total streamwise velocity in the x direction is U(y,z,t) + u(x,y, z,t), the cross- 
stream velocity in the y direction is V(y, z, t) + v(x,y, z,t), and the spanwise velocity in the z direction is 
W(y, z, i) + w(x,y, z,t). The flow is confined to the channel \y\ < 1, \z\ < L z /2 and the Reynolds number, R, 
is based on the half cross-stream channel distance L y = 1. With constant velocity channel walls forcing a Couette 
flow as a laminar equilibrium the boundary conditions are U(±l, z,t) = ±1, V(±l,z, t) = W(±l,z,t) = and 
u(x,±l, z,t) = v(x, ±1, z, t) = w(x, ±1, z, t) = which imply that the normal vorticity 77 = d z u — d x w satisfies 
the boundary condition 77(21, ±1, z, t) = and from continuity, v y (x, ±1, z, t) = 0. Periodic boundary conditions are 
imposed in the spanwise direction. 

We write perturbation equations about the streamwise mean flow, U(z,y,t), neglecting the small spanwise and 
cross-stream mean flows, as it can be verified that ||V^||,||VF|| << \\U\\. With this simplification the perturbation 
equations can be reduced, following steps similar to those used in the derivation of the Orr-Sommcrfcld and Squire 
system [17, 54] , to two equations in the normal velocity, v, and normal vorticity, 77. The resulting equations for the 
perturbations in these variables are: 

v t + A- 1 (UA X + U zz v x + 2U z v xz - U yy v x - 2U z w xy - 2U yz w x - AAv/R) = F v , (2a) 

Vt + Urj x - U z v v + U yz v + U v v z + U zz w - Ar]/R) = F n . (2b) 
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where A -1 is the inverse of the Laplacian, A = d xx + d 2 y + d zz . On the RHS F v and F v represent deviations of the 
perturbation-perturbation advection from its streamwise average, as is the convention in Reynolds averaging, together 
with input from FST. A parameterization of these terms will be specified in the next section. 
The spanwise and streamwise velocity expressed in terms of the variables v and r\ are: 

A h w = -v yz - rj x , A h u = -v yx + r) z , (3) 

in which A/j = d xx + d zz denotes the horizontal Laplacian. 
The mean streamwise flow, U , evolves according to 

U t = -{UV + uv) y - (UW + uw) z + AU/R . (4) 

The mean streamwise flow U is maintained against dissipation by the streamwise component of the force from the 
eddy Reynolds stresses and by the acceleration induced by the mean roll circulation —{UV) y — (UW) Z which can be 
written equivalcntly as —U y V — U Z W . In spanwise independent flows this term reduces to — U y V and represents the 
familiar lift up mechanism [6] . In order to evolve the streamwise mean flow the eddy Reynolds stresses and the fields 
V and W associated with the roll circulation arc required. Wc define the streak perturbation as the component of 
U that deviates from its spanwise average, [U] = Udz/L z , so that that the velocity associated with the streak is 
U s = U-[U]. 

The streamwise mean cross-stream and spanwise velocities; V and W respectively, can be obtained from the mean 
streamfunction ^(y, z, t): 

V = -^> z , W = ^ y , (5) 

which evolves according to the equation for mean streamwise vorticity A\I/. 

A* t = (VW + vw) zz - [VW + vw) yy ~ (W 2 - V 2 + w 1 - «2) y2 + AAf/i? . (6) 

This vortex is forced by the torque arising from perturbation Reynolds stresses as discussed in the previous section. 
The eddy torque is the only forcing maintaining the roll circulation against dissipation as the mean terms in (6) 
produce only advection of streamwise vorticity in the cross-stream/spanwisc plane. 

Equations (2a, 2b, 4, 6) comprise the roll/streak/turbulence dynamics. In the absence of FST this system has as a 
stable equilibrium solution only the perturbation stable Couette flow U = y, V = W = 0. In the presence of FST the 
combined roll/streak/turbulence dynamics includes the ensemble mean Reynolds stress from the perturbation field 
described by (2a, 2b, 4, 6) giving rise to new equilibria. We next show how SSST can be used to find these new stable 
equilibria. 



IV. THE SSST SYSTEM GOVERNING ROLL/STREAK/PERTURBATION DYNAMICS 

The SSST system[43] includes the three components of the streamwise mean flow (4, 6), and the ensemble mean 
Reynolds stress from the perturbations (2a, 2b). In the perturbation equations (2a, 2b) stochastic excitation is 
introduced to parameterize both the exogenous FST and the endogenous scattering by eddy-eddy interactions. The 
perturbation equations with this parameterization comprise the Stochastic Turbulence Model (STM). The STM 
provides accurate eddy structure at the energetic scales because the highly non-normal dynamics associated with the 
non- normal linear operator in strongly sheared flows predominates in determining the perturbation structure [28, 55]. 
This parameterization has been widely used to describe the dynamics of turbulence in channel flows [24, 26, 27, 29-33] 
and has also been instrumental in advancing robust control of channel flow turbulence [28, 56-59]. The STM has also 
been verified to determine with great accuracy the eddy structure of the midlatitude atmosphere [48, 49, 60-64]. 

We use Fourier expansion in the streamwise direction, x, for the perturbations that deviate from the streamwise 
mean: 

v = Y,v k (y,z,ty kx , r 1 = Y J Uv,z,t)e ikx , (7) 

k k 

in which the k = streamwise wavenumber is excluded. We discretize the perturbation equations (2a, 2b) in the cross- 
stream, y, and spanwise, z, directions with the state 4>k = [vk, Vk] T prescribed by the normal velocity and vorticity on 
a y — z grid for each x Fourier component. Streamwise mean Reynolds stresses can be obtained from the streamwise 
mean covariance matrix of the perturbation state. Under the ergodic assumption this streamwise mean covariance is 
the same as the ensemble mean covariance [89], Cfc =< <pk</>l > ( in this expression < ■ > denotes ensemble averaging 
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and the subscript indicates that the statistics are those of the eddy field components with streamwise wavenumber 
k). 

In order to evolve the perturbation covariance we must first specify an excitation to maintain the FST. We take an 
excitation in (2a, 2b) that is delta correlated in time and of the general form: 



where is a normally distributed independent random column vectors of length equal to twice the number of 
discretization points, that satisfies: 

<&(*)& 00 >=S ij 6(t-s), 

with the structure matrix F determining the spatial coherence of the forcing of the cross-stream velocity and vorticity 
i.e. the fluid is excited delta correlated in time with a spatial structure that is a superposition of the columns of F. In 
the ensemble equations for the quadratic eddy covariance at streamwise wavenumber k the stochastic forcing enters 
through its covariance Q k : 

Q fc = FF f • (8) 

The qualitative features of the SSST dynamics are insensitive to the structure of the forcing, Q^, as long as it 
excites the most energetic structures. The reason is the high non-normality of the perturbation operator which leads 
to strong amplification of a few optimal structures that in turn determine the perturbation field[24, 26]. 

We take Qfc to be proportional to M^ 1 , where is the metric that determines the perturbation kinetic energy 
at streamwise wavenumber k through the inner product: 

E k = $lM k 4> k . (9) 

This forcing covariance excites the system so that each degree of freedom receives equal energy. The energy metric is 
given by 

M fc = ± (I#L* + L*tL* + I#L*) , (10) 

where ii k = L^^fc, v k = L^fej and w k = li^(f> k . Explicitly, the linear operator Ljj is the projection, 

L* = [I 0] , (11) 
while the two other linear operators are obtained using Equation (3): 

fc _ / -ikA^d v \ k _( -A^d 2 yz o \ 
K -{ -A^d z ) ' L » - (, -ikA- h " ) • (12) 

The ensemble averaged covariance evolves according to the deterministic Lyapunov equation [47]: 



—± = A k (U)C k +C k A[(U) + f 2 Q k , (13) 

in which f 2 is an amplitude factor and A k (U) is the linear operator in (2a, 2b) which depends on the streamwise flow 
U(y,z,t). In matrix form the operator in (13) is: 



A k (U) = 




(14) 



with 

Los = A" 1 (-ikUA + ik(U yy - U zz ) - 2ikU z d z - 2ik(U z d% yz + U yz d 2 yz )A^ + A 2 /R) , (15a) 

L Cl = 2fc 2 A" x (U z y + U yz ) A~\ (15b) 

Lc 2 = U z d y - U y d z - U yz + U zz 8 2 yz A-\ (15c) 

Lsq = -ikUA + ikU zz A' h 1 + A/R (15d) 
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The covariances, C k , evolved by (13) provides the Reynolds stresses for the mean flow equations (4, 6). For example 
the Reynolds stress uv is given by: 



^Re Uag ff>£C fci L^ 



(16) 



where diag denotes the matrix diagonal and n the number of streamwise harrmonics. All the Reynolds stresses can 
be written similarly as linear functions of the covariance matrix and the streamwise mean equations (4, 6) can then 
be expressed concisely in the form: 

^ = G(r) + LC (17) 

where T = [U, ^>] T denotes the streamwise mean flow, G a function of the mean flow that includes the dissipation and 
external forcing, C = [C^, • • • , CfcJ and LC is the forcing of the mean by the Reynolds stresses, with L a linear 
operator. 

Equations [(13), (17)] comprise the SSST system for the roll/streak/turbulence dynamics: 



dC 



k 



(It 

dT 
~dt 



A fe (Pr)C fc +C fc A[(Pr) + / 2 Q fc , (18a) 
G(r) + LC, (18b) 



with P the projection of V onto the mean streamwise flow so that Pr = U. The SSST dynamics can be written 
concisely as 

^ = S b) > ( 19 ) 

by defining the SSST state X = [C, T] T . The equilibr ium states satisfy S(xeq) = 0. 

Equation (19) constitutes a closed, deterministic, autonomous, nonlinear system for the co-evolution of the stream- 
wise mean flow and its consistent field of turbulent eddies. Although the effects of the turbulent fluxes are retained in 
this system, the fluctuations associated with turbulent eddy dynamics are suppressed so that the dynamics of turbu- 
lent eddy /mean flow interaction and particularly the equilibria arising from this interaction are revealed with clarity. 
These nonlinear equilibrium states are intrinsically associated with the turbulence and are therefore dynamically dis- 
tinct from coherent nonlinear states that have been extensively studied in geophysical [66, 67] and in laboratory shear 
flows [68-70]. 

The SSST system is globally stable [43] and the attractor of the SSST system may be a fixed point, a limit cycle, 
or a chaotic attractor. Examples of each of these behaviors has been found in the SSST description of geophysical 
and plasma turbulence [41, 43, 45]. 

The concept of the SSST system trajectory is novel because it is not the trajectory of a realization of the turbulent 
system but rather the trajectory of the statistical mean state of the turbulence which evolves on the time scale of the 
mean flow. 

The SSST system introduces a new stability concept to fluid dynamics which is the stability of an equilibrium 
between a mean flow and its consistent field of turbulence. This stability theory generalizes the traditional hydro- 
dynamic stability theory of Rayleigh[71]. If a mean flow is perturbation unstable (in the sense of Rayleigh) it is 
also structurally unstable (in the sense of SSST). However, the converse is not true and perturbation stability does 
not imply structural stability. In fact, emergence of roll/streak structures in shear flow will be shown to occur in 
association with structural instability of a perturbation stable state. 



V. STABILITY ANALYSIS OF ROLL/ STREAK/ TURBULENCE EQUILIBRIA 



Assume that for the given forcing covariance, / 2 Qfc, the equilibrium \eq = [C eg ,r e(? ] T of the SSST equations 
(18) has been determined. We can study its stability by linearizing the SSST system about this equilibrium. The 
perturbation equations take the form 



dSC k 
dt 
dST 
~~dT 



A k (PT eq )8C k + SC k A{(Pr eq ) + 5A k C keq 



C keq 5A{ 



k ' 



dG 
dT 



sr + lsc 



(20a) 
(20b) 
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where ST = [5U, 5ty] T is the perturbation in the mean flow quantities: the streamwise flow, 6U, and the roll stream- 
function, Sty. The perturbation to the operator, that controls the eddy field is SA.^. This operator perturbation 
is produced by perturbation to the mean streamwise flow, SU. Setting 8\ = [SC, ST] T the perturbation equations can 
be written concisely as: 

The linear operator C = dS/dx\ Xcq depends on the equilibrium state \eq = [C eq ,T eq ] T which in turn depends on the 
Reynolds number R, the channel geometry, the mean flow forcing, and the stochastic excitation / 2 Qfc- Eigenanalysis 
of the linear operator C then determines the structural stability of this equilibrium roll/streak/turbulence complex [90]. 

The familiar laminar Couette flow equilibrium T eq = [U eq = y, ty eq = 0] T and C eq = is a solution of the SSST 
equations (18) in the absence of FST (/ = 0). Because at zero forcing the eddy covariance vanishes, C eq = 0, 
the first of the perturbation SSST equations (20) reduces to an unforced Lyapunov equation which inherits the 
perturbation stability of Couette flow at all Reynolds numbers [73], i.e. the stability of A(Pr eg ). The second equation 
is asymptotically unforced and clearly stable so in the absence of FST the system (20) is structurally stable as well 
as perturbation stable. From this argument it is clear that structural instability of a flow that is perturbation stable 
in the sense of Rayleigh requires non- vanishing C eq or equivalently non zero values of /. 

In the presence of spanwise homogeneous FST there is a class of spanwise independent equilibria T eq = 
[U eq (y),ty eq = 0] T with non vanishing C eq . In these equilibria the mean streamwise flow U eq (y) is maintained 
by a balance between diffusion and the component of Reynolds stress divergence —(uv) y in the inhomogeneous (cross- 
stream) direction. The equilibria are possible because all the Reynolds stresses are independent of z, and symmetry 
requires that WW = 0. These spanwise independent ensemble equilibria in the presence of FST correspond to boundary 
layer flow equilibria that depart from the Couette profile in y but have no z dependence. We will demonstrate that 
for sufficient amplitude of FST these equilibria, while remaining perturbation stable, become structurally unstable 
giving rise to roll circulations with associated streaks. 

For convenience the amplitude of the forcing, is chosen to maintain RMS perturbation velocity 1% of the mean 
Couette flow velocity when it is used to excite the Couette flow so that when this excitation is introduced into the 
perturbation variance equation as f 2 Qk the adjustable amplitude, /, corresponds approximately to RMS FST as 
a percentage of the mean flow velocity. RMS perturbation velocity is very nearly linearly proportional to / as / 
increases prior to the bifurcation to roll /streak equilibria although deviating slightly because the mean flow profile 
deviates from Couette as FST increases (cf. Fig. 6). 

We demonstrate the structural instability of the spanwise independent equilibria T eq = [U eq (y), ty eq = 0] T in the 
presence of FST concentrated in a single wavenumber fc. This simple case reveals the character of the instability 
and the results do not change qualitatively when multiple wavenumbers of FST are included. We first examine an 
example at R = 400, wavenumber k = 1, and spanwise periodic channels on the interval 3/20 < L z /(2ir) < 3/10. 
These channels are narrower than the minimal channel of Hamilton, Kim and Waleffe[74] (heraftcr HKW) for which 
L z /(2tt) = 6/10. The calculations were performed with N y = 21 and N z = 20 points. Convergence was verified by 
repeating the calculations at higher resolutions. 

We use the power method to find the structure and growth rate of the most unstable eigenmode of the C operator 
in (21) for the spanwise independent equilibrium flow T eq = [U eq (y), ty eq = 0] T . Contours of the growth rate of the 
most unstable eigenmode as a function of STM excitation parameter / and channel width L z are shown in Fig. 2. The 
maximum growth rate of the C operator increases with / and for a critical intensity, f c (L z ), the spanwise independent 
flow becomes structurally unstable with eigenmodes in the form of exponentially growing roll/streak structures. 
The growth rate decreases with channel width at constant / and for sufficiently narrow channels L z /(2ir) < 0.205 
the spanwise independent mean flow is structurally stable for all / and no roll/streak equilibrium is supported, in 
qualitative agreement with the findings of the minimal channel simulations of Jimenez and Moin[75]. 

For a channel with L z /(2tt) = 0.3 the spanwise independent equilibrium is structurally unstable for / > f c = 8.25. 
The growth rate of the most unstable eigenfunction for this channel with STM excitation parameter / = 12.8 is 
A = 0.0166. The structure of the most unstable eigenfunction is shown in Fig. 3. The eigenfunction comprises both 
a mean flow perturbation, ST, and an eddy covariance perturbation, SC. The Reynolds stresses associated with SC 
produce accelerations and torques in exact agreement with the mean flow perturbation consistent with exponential 
growth. These structural instabilities typically equilibrate to finite amplitude roll/streak equilibria similar in structure 
to the most unstable eigenfunction, as seen in Fig. 4. The critical f c at which structural instability occurs is a 
bifurcation point in a diagram of equilibria as a function of /, as shown in Fig. 5, 6. 

We denote by E r the mean kinetic energy of the roll, obtained by averaging (V 2 + W 2 )/2 over the channel; by Ejj 
the mean kinetic energy of the streamwise mean flow, obtained by averaging U 2 /2; by E s the mean kinetic energy 
of the streak, obtained by averaging U 2 /2; and by E p the mean kinetic energy of the eddies, obtained by averaging 
(u 2 + v 2 + w 2 )/2. As / increases the STM excitation parameter rises and the increasing Reynolds stress induces 
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departure of the cross-stream flow from the pure Couette flow but for / < 8.25 the mean flow remains uniform in 
the spanwise direction and these spanwise independent equilibria are stable fixed points of the SSST system. The 
RMS velocity of the perturbation field y/E p /2 is shown as a function of / in Fig. 6. Note that as the spanwise 
uniform equilibria bifurcate to roll/streak equilibria for / > f c = 8.25 the rate of increase of perturbation RMS 
velocity decreases as the turbulence is diverted to drive the roll/streak structure which, being lightly damped, reaches 
high equilibrium velocity. It is interesting that the RMS velocity remains approximately 10% of the background flow 
velocity until forcing excitation f u = 13.5 at which point the roll/streak equilibrium undergoes a second structural 
instability. The roll/streak equilibria for f c < f < f u arc perturbation stable, and the breakdown that occurs at f u is 
due to a secondary structural instability of the finite amplitude roll/streak. This secondary structural instability will 
be examined in future work. 

We next demonstrate structural instability of the Couette flow in the minimal channel considered by HKW[74] 
taking L z /(2ir) = 0.6, R = 400 and the smallest streamwise wavenumber in the HKW channel, k = 1.143. These 
calculations were made with N y = 21 and N z = 40 points. This flow bifurcates from spanwise independent equilibria 
to spanwise dependent equilibria at f c = 5.82 as shown in Fig. 7. Note the qualitative similarity with the bifurcation 
diagram for the smaller channel shown above. Both the RMS streak velocity and the RMS roll velocity vary as 
V/ — fc near the the bifurcation point (dashed curve in Fig. 7) consistent with a supercritical pitchfork. At f u = 8.45 
the roll-streak equilibrium loses structural stability and no nearby equilibrium or periodic solution exists for / > f u . 
For values close to this second structural instability the equilibria exhibit a \J f u — f behavior consistent with a second 
order subcritical bifurcation which will be examined in future work. The most unstable perturbation of the SSST 
system about the unstable equilibrium state without roll/streak structure is shown for STM excitation parameter 
/ = 8.4 < f u in Fig. 8. 



VI. STRUCTURE OF THE ROLL/STREAK EQUILIBRIA 

The streamwise mean flows of the structurally stable roll/streak equilibria in the HKW channel for STM excitation 
parameter in the interval f u > f > f c are shown in Fig. 9. These equilibria exhibit streamwise high and low speed 
streaks that increase in amplitude as / increases. A useful measure of streak strength is its lift angle [23] which we 
define here as: 

_! / d z U \ 
s =maxtan- , 

where [U] is the spanwise average streamwise mean flow U. These equilibrium streaks reach 6 S = 56° at / = 8.45. 
Despite their high lift angles all these equilibria are perturbation stable. The eigenvalues, a, of the perturbation 
operators, A&, for the flows in Fig. 9 are shown in Fig. 10. Note the emergence of a mode with frequency Oj = as 
the streak increases in amplitude with increase in /. This is the sinuous mode that is associated with the spanwise 
inflection of the streak and that is often assumed to be responsible for streak breakdown[54, 74, 76-78] . However at 
/ = 8.45 the roll/streak flow is still robustly stable and the instability that occurs at f — 8.45 is solely a structural 
instability of the cooperative turbulence/mean flow SSST dynamics. 

Despite the large lift angle the equilibrium streamwise mean flows do not resemble the mean flows of the turbulent 
state. In Fig. 11 we compare the spanwise averaged streamwise mean flow [U] for the flows with STM excitation 
parameter / = (8.4, 7.5, 6) with the corresponding time averaged mean flow in the turbulent HKW channel. The 
equilibrium mean flows that occur in this range of STM excitation parameter indicate that these roll/streak equilibria 
are laminar. Indicative of this laminar regime is the viscous dissipation of the streamwise mean flow: 

D = ^J Q L dz fl d y( u v + u * + v v + ^ 2 + w v + w *) > 

The ratio D/Dc, where Dc is the dissipation associated with the Couette flow, of the equilibria for FST intensities 
8.45 > / > f c = 5.82 is in the range 1 < D/Dq < 1.4 while this ratio is of order 3 in the turbulent state. 

The laminar roll/streak equilibria shown in Fig. 9 have spanwise wavenumber 2. The spanwise width of the channel 
is 90y + for the equilibrium with / = 8.4 (the wall unit is defined as y + = 1/ y/R\U y ], where R is the Reynolds number 
and [U y ] is the mean shear at the boundary) implying streak spacing 45y + which is half that found in turbulent 
boundary layers. However, it should be kept in mind that this wall unit is being calculated for an essentially laminar 
flow. In section VIII we show that this spacing does agree with the observed streak spacing of 2 displacement 
thicknesses that is observed in laminar boundary layers before transition [79]. 

Consider the mechanism maintaining the roll/streak structure at STM excitation parameter / = 8.4 (Fig. 9d). The 
roll circulation is maintained against friction only by the turbulent stress divergence in (6) as the quadratic streamwise 
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mean terms in (6) do not generate mean streamwise vorticity. The cross-stream/spanwisc acceleration (V, W) due to 
the eddy flux divergence is shown in Fig. 12c. Note that this acceleration is consistent with the circulation shown 
in Fig. 13a. The acceleration induced by the mean momentum flux divergence is subdominant in these equilibrium 
solutions and as a result the total acceleration (V, W) has the structure of the acceleration induced by the eddies as 
shown in Fig 13c and Fig 13a. 

While the roll circulation is maintained against friction solely by the torque induced by the Reynolds stress diver- 
gences, the streak is influenced both by Reynolds stress divergence and by mean momentum flux divergence (cf. Eq. 
4). The mean momentum flux divergence can be identified with the lift up mechanism as shown in Fig. 12b and this 
mechanism dominates in the streak maintenance (cf. Fig. 12b and Fig. 12d). The eddy Reynolds stress divergence, 
shown in Fig. 12d, tends to damp the streak consistent with the eddies extracting energy from the spanwise shear. 
This mean deceleration of the streak by the 3-D eddies in laboratory shear flow contrasts with the acceleration by 
quasi 2-D eddies that is primarily responsible for jet formation in planetary atmospheres [80-82]. In compensation for 
the loss of this dominant 2-D jet formation mechanism, these 3-D shear flows gain an indirect pathway for maintenance 
of the streak: the Reynolds stress divergences induce roll circulations which through the lift-up mechanism maintain 
the streaks against both viscous dissipation and the deceleration induced by the turbulent Reynolds stress divergence. 
This dual role of the eddy field in maintaining the equilibrium roll/streak structure will be discussed further in the 
next section. 

We turn now to the structure of the eddy field at equilibrium. Eddy structures can be ordered in energy by 
cigenanalysis of M^CM 1 ' 2 , where M 1 / 2 denotes the square root of the metric M. Eigenfunctions of M 1 / 2 CM 1 / 2 
in descending order of eigenvalue define the empirical orthogonal function (EOF) or Karhunen-Loeve decomposition 
of the eddy field. The percentage of the energy accounted for by each of the first 40 gravest modes of the covariance is 
shown in Fig. 14 for STM excitation parameter / = 5, for which the flow is spanwise uniform, and also for intensity 
parameters / = 6, 7.5, 8.4. It is clear that the variance is spread over many structures but as / increases the first 
EOF begins to dominate and its structure becomes a good representation of the eddy structure. This dominant EOF 
for the equilibrium structure at / = 8.4 is shown in Fig. 15. This eddy field is characterized by sinuous oblique waves 
centered at the wings of the streak and slanted in the vertical. This is the structure that produces the coherent torques 
maintaining the roll circulation. The dominant EOF is close in structure to the least stable mode of the flow which is 
shown in Fig. 16. Note however that this mode is robustly stable so it is maintained by excitation and non-normality, 
with the latter dominant. The large excitation of the mode is due to the interaction between this mode and the other 
modes of the system as revealed by its optimal excitation structure which is its adjoint [7, 47, 83]. The adjoint of the 
least stable mode in the energy inner product is shown in Fig. 17. An initial condition consisting of its adjoint excites 
the least stable mode at amplitude a factor 1900 greater than an initial condition consisting of the least stable mode 
itself so the mode arises out of the FST primarily due to excitation of its adjoint. 

Because the excitation is chosen to be white in energy and all modes are stable the structure of the eddy field can 
be understood dynamically by examining optimal structure evolution. The optimal perturbation that leads to the 
greatest growth in energy in 10 time units for the equilibrium at / = 8.4 is shown in Fig. 18. The energy growth 
of this optimal is close to the energy growth of the adjoint of the least stable mode (cf. Fig. 19). Evolution of the 
maximum mean streamwise torque, G x , induced by the Reynolds stress divergence of these optimal perturbations is 
also shown in Fig. 19. The torque increases as the perturbation energy increases. The energy evolution of the t = 10 
optimal for the equilibria with / = 7.5 and / = 5, for which value the flow is spanwise independent, are also shown. 
The structure of the optimal when it reaches its maximum energy at t = 15 is shown in Fig. 20. The structure of this 
evolved optimal perturbation is similar to the structure of evolved optimals in equilibria at lower FST intensities as 
expected from the universality of the dynamics of oblique perturbations in three dimensional shear flows [1, 13, 14, 84]. 
The spanwise streak serves to collocate the perturbation structures aligning them so that the spanwise Reynolds stress 
divergence produce torques in phase with the evolving roll. 

Finally we note that when the total field of the equilibrium mean flow and a typical realization of the eddy field are 
plotted together only a weak undulation of the streak structure can be discerned as is observed with laminar streaks 
before transition. For example a sample realization of the total equilibrium flow at STM excitation parameter / = 8.4 
is shown in Fig. 21. 

VII. MECHANISM OF ROLL/STREAK EQUILIBRATION 

We wish to gain an understanding of the dynamics underlying equilibration of the structural instability of the 
spanwise independent flow in part to provide insight into turbulent equilibria in general and in part as a first step in 
understanding how loss of structural stability by these roll/streak equilibria at / = 8.45 leads to transition to to a 
time dependent state. 

We will describe the equilibration of the structural instability of the spanwise independent flow at STM excitation 
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parameter / = 8.4 and in particular show how the inflectional mode is instrumental in producing the equilibrium. 
When the most unstable SSST eigenfunction is introduced into the spanwise independent flow, it grows exponentially 
at rate A = 0.014 as predicted by SSST theory and finally equilibrates as can be seen in Fig. 22. As the STM 
excitation parameter approaches / = 8.4 oblique non-normal wave perturbations dominate the structures forcing 
the roll through their Reynolds stress and these produce a roll/streak equilibrium with highly inflected streaks (cf 
Fig. 9d). In association with this inflection the primary inflectional mode (with structure as in Fig. 16) approaches 
the stability boundary (cf Fig. lOd). Because it is drawing energy from the streak this mode produces strong 
downgradient Reynolds stresses that damp the mode while the Reynolds stress it produces to force the roll circulation 
are relatively weak. In order to compare the relative contribution of the direct downgradient Reynolds stress due 
to the inflectional mode in damping the streak with the mode's indirect Reynolds stress effect in forcing the roll 
circulation and thereby building the streak we artificially impose a modification of the real part of the eigenvalue 
of the mode at equilibrium, a r E, specifically we set the growth of this mode equal to 9/lOay^ and 11/lOoy^;, and 
integrated forward the SSST equations in order to determine the mean flow tendency. When the mode is less stable 
the perturbation energy increases, the associated roll circulation also increases as it is directly forced by the oblique 
structure of the perturbations, but the streak decreases because the enhanced downgradient fluxes by the less damped 
mode dominate over the increase in the streak induced by the roll circulation. The opposite happens when the mode 
is made more stable, as shown in Fig. 23. 

We conclude that the primary mechanism of streak stabilization at high FST is the inflectional mode. As the / 
rises above / = 8.45 the inflectional mode is no longer able to stabilize the streak and a second structural instability 
ensues at f u = 8.45 in which the oblique waves further accelerate the roll/streak complex. However, the streak 
remains perturbation stable until a very high amplitude is reached at which point the flow becomes time dependent 
and aperiodic so that the notion of an unstable temporal mode is no longer well defined. We will not examine this 
time dependent regime further in this work. 

VIII. MECHANISM PRODUCING THE SPANWISE STREAK SPACING 

A fully turbulent boundary layer, such as that approximated by the Reynolds-Tiederman profile, is maintained by 
the ensemble of eddies in the boundary layer. It is commonly observed in turbulent boundary layers that spanwise 
streak spacing is approximately 100 y + with wall unit y + = y v] (dU/dy) and dU / dy being evaluated at the boundary. 
As the boundary layer is itself approximately 50 y + in wall normal extent this spacing is consistent with a roll of unit 
aspect ratio confined to the boundary layer[85, 86]. 

We have concentrated on the formation of streaks from FST in which the deviation of the mean boundary layer 
flow profile from the stationary Couette flow is small compared with that found in fully turbulent boundary layers. 
In order to study streak spacing in a numerically resolved example we maintain a Blasius profile stationary with an 
appropriate body force [79], subject it to a supercritical STM excitation parameter, and obtain the maximum growth 
rate of the structural instability as a function of spanwise wavenumber m of the unstable streak using the power 
method. Growth rate of the most unstable SSST roll/streak eigenmodes at STM excitation parameter / = 10 in 
a Blasius boundary layer at Reynolds number R x = 1.6 x 10 5 (based on the distance from the leading edge) are 
shown in Fig. 24. The maximum SSST instability occurs at the wavenumber corresponding to a spacing between low 
speed streaks of Az = 2.4#i where Si = \.11\J vxjU is the displacement thickness, consistent with unit roll aspect 
ratio. Although the implied selectivity is not very strong this result also agrees with observations [79]. The maximum 
is achieved by oblique waves with obliqueness parameter = tan -1 m/k, close to the value 53° obtained from the 
simple argument of section I. This agreement shows that the basic dynamics are captured by the oblique plane wave 
solutions on an unbounded constant shear flowfl, 13, 14, 84]. These structures are scale independent and the streaks 
that are formed by the SSST instability share the universal character of these oblique perturbations. 

IX. DISCUSSION 

There are a number of points we wish to emphasize in connection with the above results: 

1. Streaks can arise from a spontaneous cooperative exponential instability of the roll/streak/shear plus turbulence 
complex. 

2. At finite amplitude the streaks of the roll/streak/shear plus turbulence instability complex typically form stable 
nonlinear equilibria. 
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3. The cross-stream/spanwisc torques supporting the equilibrium roll/streak structure are produced by oblique 
waves properly collocated with the roll/streak structure. 

4. The structural stability boundary for streak formation by this cooperative instability is not associated with 
modal instability of the streak and the finite amplitude equilibrium roll/streak structure is stable despite being 
highly inflectional in the spanwisc direction. In fact, the inflectional mode acts to stabilize the streak. 

A question often raised is whether the mechanism forming the roll/streak complex is essentially linear non- normal 
or essentially nonlinear [54, 87]. We find that both are involved in agreement with Kim and Lim[22]: the roll/streak is 
the optimal non- normal growth structure the growth of which is essentially related to non- normality of the linear shear 
flow dynamics while the dcstabilization of this roll/streak structure arises from the quadratically nonlinear Reynolds 
stresses. This forcing by perturbation Reynolds stresses of the non-normal optimal roll structure is an example of 
the mechanism of structure maintenance in turbulence by nonlinear scattering of perturbation energy into the linear 
non-normal optimally growing subspace. Because the non-normal roll/streak structure is so highly amplified, the 
energy scattered by nonlinear interaction into this structure dominates the structure of boundary layer turbulence. 

Hamilton et al.[74] (HKW) obtain in an integration of a spatially constrained shear flow model what they identify 
as a regeneration mechanism for maintaining at finite amplitude the roll/streak complex. Their mechanism consists 
of a dynamic process of continual streak growth and decay associated [88] with forcing of the roll circulation by an 
unstable mode arising from spanwisc inflectional streak instability. Their work addresses the problem of streak self- 
maintenance and presupposes that a streak of finite amplitude sufficient to produce an inflectional instability already 
exists. This work addresses the formation of streaks from arbitrarily small initial perturbations in FST. The finite 
amplitude equilibria we find, while inflectional, are modally stable and our streaks are maintained by the wave/mean 
flow interaction arising from a large subset of non-normal transient perturbation structures rather than by a single 
inflectional mode. In fact, as we have seen, the inflectional mode is primarily responsible for stabilizing the SSST 
streak instability at finite amplitude. Nevertheless, their numerical simulations agree with both observations and SSST 
in a number of other particulars including the importance of the strcamwise roll as the linear non-normal optimally 
growing structure necessary for the formation of streaks and of Reynolds stresses in forcing the strcamwise roll. HKW 
find it remarkable that turbulent perturbations "..produce additional strcamwise vorticity in exactly the right places 
to augment the strcamwise vortices." [74] but as we have seen, at least in the initial stages of streak formation, this 
coincidence is a necessary consequence of the existence of the SSST eigenmode and its nonlinear extension. 

The SSST system provides analysis tools for obtaining a fundamental understanding of the underlying mechanisms of 
turbulence that can not be gained from interpretation of simulations alone so that we may say that SSST constitutes a 
theory of turbulence as distinct from a description of turbulence or an interpretation in dynamical terms of observations 
of turbulence. In this work we have used SSST to study the formation and maintenance of roll/streak structures in 
FST which is important in its own right but also as a component of the mechanism of bypass transition which typically 
proceeds from a pre-existing roll/streak structure. 

X. CONCLUSION 

Emergence of the roll/streak coherent structure in turbulent flow is a problem of great theoretical and practical 
importance. In this work we applied SSST to demonstrate a mechanism by which the roll/streak structure can arise 
from an emergent instability of roll/strcak/turbulence interaction in boundary layer shear flow. 

This emergent SSST instability giving rise to the roll/streak structure exploits the optimality of the non- normal 
roll/streak structure growth mechanism not by introducing an individual chance perturbation in cross-stream velocity 
but rather by organizing the ubiquitous torques associated with turbulent Reynolds stress divergence in the cross 
stream/spanwise plane to produce and maintain the optimal roll structure. This organization of the Reynolds stress 
by the streak resulting in forcing of the roll provides the missing coupling between the streak and roll that is required 
to produce instability from the non- normal growth process. At small amplitude the roll/streak structure grows 
exponentially as an eigenmode but at finite amplitude this growth is arrested and the structure approaches a nonlinear 
equilibrium. This mechanism of streak formation and equilibration is not related to instability of the perturbation 
dynamics of the streak and this cooperative SSST instability occurs in the absence of shear flow instability and the 
inflectional mode is found to stabilize the streak. 
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FIG. 1: (Color online) Left: Streamwise flow perturbation 5U = ecos(jiy / {2L y )) cos(2irz/Lz) imposed on a background 
Couette flow to examine the mechanism of turbulent Reynolds stress organization. Right: Resulting equilibrium Reynolds stress 
divergence induced tendency in the cross-stream/spanwise streamfunction, dty/dt, normalized by the mean flow perturbation 
amplitude, e, and the square of the STM excitation parameter parameter, /. Imposition of a spanwise perturbation breaks the 
spanwise symmetry of the Couette flow producing a coherent streamwise torque proportional to both the mean flow perturbation 
and to the eddy field variance. The channel dimensions are L y = 1, L z = 1.2n, the Reynolds number is R = 400, and the eddy 
field is at streamwise wavenumber k = 27r/(1.757r). 
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FIG. 2: (Color online) Growth rate of the most unstable eigenfunction of the SSST system linearized about the spanwise uniform 
equilibrium as a function of STM excitation parameter, /, and spanwise channel width, L z . Channel width L z /(2n) = 3/10, 
as used in the example of Fig. 3, lie on the abscissa of this plot. Channels with spanwise width L z /(2tt) < 0.205 are stable. 
The perturbation field comprises a single wavenumber, k = 1. The Reynolds number is R — 400. 
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FIG. 3: (Color online) The most unstable eigenfunction of the SSST system linearized about the structurally unstable equilib- 
rium with no spanwise variation at STM excitation parameter / = 12.8. The growth rate of this eigenfunction is A = 0.0166. 
(a): streamwise mean cross-stream/spanwise velocity vectors (SV,SW) in the cross-stream/spanwise plane, (b): streamwise 
mean streamwise velocity 8U associated with the same eigenfunction (negative values dashed). The ratio of the maxima of the 
fields (SU, SV, SW) is (1,0.06,0.03). The unstable eigenfunction also has a perturbation covariance component, SC, the effect 
of which is indicated by the acceleration these perturbations induce in the corresponding velocities (cf. equations (4) and (6)): 
(c): (5V,5W). (d): 5U. Parameters are k = 1, Lz/(2ir) = 0.3 and R = 400. 
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FIG. 4: (Color online) Finite amplitude equilibrium roll/streak structure at STM excitation parameter / = 12.8. Shown are 
contours of the streak velocity, U s , and velocity vectors (V, W) of the roll circulation that equilibrates from the most unstable 
eigenfunctions shown in Fig. 3. 
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FIG. 5: (Color online) Roll/streak equilibria as a function of the STM excitation parameter, /. At the critical turbulence level 
corresponding to f c = 8.25 the spanwise independent equilibria bifurcate to spanwise dependent roll/streak equilibria. Shown 
are RMS streak amplitude (continuous), RMS roll amplitude (dashed), and RMS streamwise mean flow amplitude (dash-dot) 
normalized by the RMS velocity of the unperturbed Couette flow, \/Ec- The Reynolds number is R = 400 and the perturbation 
field wavenumber is k = 1. 
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FIG. 6: (Color online) RMS velocity of the sum of perturbation, roll, and streak components \J (E v + E m + E s )/ Ec (dashed), 
and of the perturbation alone ^jE v jEc (solid) normalized by the Couette RMS velocity, y/Ec, as a function of STM excitation 
parameter, /. These curves diverge at the bifurcation STM excitation parameter, f c = 8.25, at which the roll and the streak 
emerge. As the STM excitation parameter increases beyond / c the roll/streak complex rapidly increases in contribution to the 
total energy while the perturbation RMS velocity remains near 10% of the mean velocity. Parameters as in Fig. 3. 
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FIG. 7: (Color online) Bifurcation diagram of the roll/streak structure as a function of STM excitation parameter, /. At 
f c = 5.825 the spanwise uniform equilibrium bifurcates to a spanwise dependent roll/streak equilibrium. Shown as a function 
of / is normalized streak strength, WO^Es/Ec) 1 ^ 2 (squares) and normalized roll strength, 10 3 (E r / Ec) 1//2 (circles). The dashed 
line indicates the y/ f — f c dependence of a second order bifurcation. The stable equilibria extend up to f u = 8.45; at which 
point the perturbation stable roll/streak equilibrium becomes structurally unstable. The dashed line indicates the y/f u — f 
dependence of a second order bifurcation. Parameters correspond to the HKW channel: length L x — 1.75-7T, spanwise width 
L z = 1.2-7T, half cross-stream height L y — 1.0 and the Reynolds number is R = 400. The perturbation streamwise wavenumber, 
k = 1.143, corresponds to the gravest mode in the channel. 
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FIG. 8: (Color online) The most unstable eigenfunction of the SSST system linearized about the structurally unstable equi- 
librium with no spanwise variation at STM excitation parameter / = 8.4. The growth rate of this mode is A = 0.014. Left: 
streamwise mean cross-stream/spanwise velocity vectors (SV, SW) in the cross-stream/spanwise plane. Right: streamwise mean 
streamwise velocity SU. The ratio of the maxima of the fields (SU, SV, SW) is (1,0.06,0.03). Parameters correspond to the 
HKW channel: length L x = 1.75-7T, spanwise width L z = 1.2-7T, half cross-stream height L y = 1.0 and the Reynolds number is 
R = 400. The perturbation streamwise wavenumber, k = 1.143, corresponds to the gravest mode in the channel. 
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FIG. 9: (Color online) Streamwise mean velocity in the (y, z) plane for the equilibria at various STM excitation parameter 
values, /. (a): at STM excitation parameter value / = 5 the equilibrium is spanwise uniform and there is no associated 
roll/streak, (b): At / = 6 the spanwise independent flow is structurally unstable and the associated equilibrium flow has a 
weak streak with associated roll/streak, (c): The equilibrium at / = 7.5. (d): The equilibrium at / = 8.4. All these equilibria 
are perturbation stable. Parameters correspond to the HKW channel: length L x — 1.75n, spanwise width L z = f .2ir, half 
cross-stream height L y — 1.0 and the Reynolds number is R = 400. The perturbation streamwise wavenumber, k = 1.143, 
corresponds to the gravest mode in the channel. 
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FIG. 10: (Color online) The least stable eigenvalues (ay, ai) of the operators that govern the perturbation stability of the 
equilibrium flows shown in Fig. 9. All the flows are perturbation stable. Note the emergence of a mode with ay = as the 
streak increases in magnitude. Parameters correspond to the HKW channel: length L x = 1.757T, spanwise width L z = 1.2n, 
half cross-stream height L y = 1.0 and the Reynolds number is R = 400. The perturbation streamwise wavenumber, k = 1.143, 
corresponds to the gravest mode in the channel. 
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FIG. 11: (Color online) The spanwise averaged streamwise flow, [U](y), for the SSST equilibria for STM excitation parameter 
/ = 8.4 (solid), / = 7.5 (dashed), / = 6 (dotted) as shown in Fig. 9. Also shown for comparison is the time and spanwise mean 
under turbulent conditions (dash-dot) as well as the laminar Couette flow. The laminar roll/streak equilibrium at / = 8.4 
produces dissipation lADc, where D c is the dissipation in Couette flow, and the half channel width in wall units at this value of 
/ is L y = 24j/ + . Parameters correspond to the HKW channel: length L x = 1.75n, spanwise width L. = 1.2-7T, half cross-stream 
height L y = 1.0 and the Reynolds number is R = 400. The perturbation streamwise wavenumber, k = 1.143, corresponds to 
the gravest mode in the channel. 
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FIG. 12: (Color online) For the equilibrium shown in Fig. 9d; (a): acceleration vectors, (V, W), of the streamwise mean roll 
circulation induced by the mean velocity field, the maximum V is 10~ 4 . (b): acceleration of the mean streamwise flow, U, 
induced by the streamwise mean roll circulation, the maximum U is 10 -2 , mainly due to the lift-up mechanism, (c): acceleration 
vectors, (V, W), of the streamwise mean roll circulation induced by the eddy field, the maximum V is 10 -3 . (d): acceleration 



of the mean streamwise flow, U, induced by the eddy field, the maximum U is 1CP' 
Parameters correspond to the HKW channel: length L x = 1.75-7T, spanwise width L z 
and the Reynolds number is R = 400. The perturbation streamwise wavenumber, k 
in the channel. 



The eddy field decelerates the streaks. 
: 1.2-7T, half cross-stream height L y = 1.0 
1.143, corresponds to the gravest mode 
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FIG. 13: (Color online) For the equilibrium shown in Fig. 9d; (a): Roll vector velocity (V, W) at equilibrium, the maximum 

V velocity is 0.02 and the maximum W velocity is 0.009. (b): Streak velocity U s = U — [U] at equilibrium, the maximum 
velocity is 0.26. (c): mean roll acceleration vectors induced by both the mean and eddy fields given by 

V x (A~ 1 ((VW+ < vw >) zz - (VW+ < vw >) yy - (W 2 + < w 2 >) yz + (V 2 + < v 2 >) V z)i). The roll circulation is 
maintained against Iriction by the eddy field, (d): mean streamwise acceleration induced by both the mean and eddy fields 
-{UV+ < uv >) y - (UW+ < uw >) z . 

The mean streamwise flow is maintained against friction by the lift up mechanism. 
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FIG. 14: (Color online) Percentage contribution of the leading EOF's of the equilibrium covariance to the total eddy mean 
energy maintained by the equilibria shown in Fig. 9. 
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FIG. 15: (Color online) Velocity field of the gravest EOF accounting for 24% of the eddy energy for STM excitation parameter 
/ = 8.4. Velocity vectors are shown with contours of streamwise velocity. Parameters correspond to the HKW channel: 
length L x — 1.757T, spanwise width L z = 1.2n, half cross-stream height L y — 1.0 and the Reynolds number is R — 400. The 
perturbation streamwise wavenumber, k = 1.143, corresponds to the gravest mode in the channel. 
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FIG. 17: (Color online) Velocity field of the adjoint in the energy inner product of the least stable mode with eigenvalue 
(ay = —0.017, tJi =0), for the equilibrium at STM excitation parameter / = 8.4. Velocity vectors are shown with contours 
of streamwise velocity. The adjoint is the optimal excitation of the mode. An initial condition consisting of the adjoint with 
unit energy excites the least stable mode a factor of 1900 greater than an initial condition consisting of the least stable mode 
itself with unit energy demonstrating that the mode amplitude derives almost entirely from non-normal growth processes. 
Parameters are as in Fig. 15. 
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FIG. 18: (Color online) Velocity field at time t — of the optimal perturbation that maximizes energy growth at t — 10 for 
the equilibrium flow with STM excitation parameter / = 8.4. Velocity vectors are shown with contours of streamwise velocity. 
Parameters are as in Fig. 15. 
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FIG. 19: (Color online) Left : Energy of the optimal perturbation that maximizes energy growth at t — 10 as a function of time 
for the equilibrium flow with STM excitation parameter / = 8.4 (solid), / = 7.5 (dashed) and for / = 5 for which there is no 
roll/streak; (dotted). Also shown is the energy growth associated with a unit energy initial condition consisting of the adjoint 
in the energy metric of the least damped mode for the equilibrium with / = 8.4 (dash-dot). Right: The time development 
of the maximum mean streamwise torque induced by the Reynold's stresses of the corresponding evolving optimals and the 
adjoint. This figure demonstrates that both the mode amplitude and its contribution to the streamwise mean torque are due 
to non-normal growth processes. All modes in these flows are exponentially stable (cf. Fig. 10a,c,d). Parameters are as in Fig. 
15. 
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FIG. 20: (Color online) Velocity field at time t — 15 of the optimal perturbation that maximizes energy growth at t = 10 for 
the equilibrium flow with STM excitation parameter / = 8.4. Velocity vectors are shown with contours of streamwise velocity. 
Parameters are as in Fig. 15. 
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FIG. 21: (Color online) Velocity field (mean and perturbation) for the equilibrium flow with STM excitation parameter 
/ = 8.4. Velocity vectors are shown with contours of streamwise velocity. The eddy field is a sum of the EOF's with amplitudes 
commensurate to their contribution to the total maintained eddy energy. The eddy field is dominated by the mean flow and 
the meandering of the streak in this laminar state is very slight. Parameters are as in Fig. 15. 
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FIG. 22: (Color online) Left: Evolution of percent RMS eddy kinetic energy 100^/ E p /Ec normalized by the RMS kinetic 
energy of the Couette flow after the structurally unstable spanwise independent equilibrium is perturbed by the most unstable 
streak perturbation shown in Fig. 8. The flow equilibrates to the roll/streak equilibrium shown in Fig. 9c and Fig. 13. Right: 
Evolution of the RMS streak energy 100 ^JeTJEc and RMS roll energy 1Q 2 ^/ E r / Ec is exponential with growth rate A = 0.014. 
Parameters are as in Fig. 8, the STM excitation parameter is / = 8.4.. 
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FIG. 23: (Color online) Percent RMS normalized streak amplitude 100 y E s /Ec (top panel), permil RMS normalized roll 
amplitude 10 3 \J E r /Ec (middle panel), and percent RMS normalized perturbation velocity amplitude 100 ^/E p /Ec (bottom 
panel) as a function of time for decrease (solid) and increase (dashed) in the damping rate of the inflectional mode compared 
to its damping rate at equilibrium. The corresponding energies at equilibrium are also shown (dotted). The inflectional mode 
clearly damps rather than drives the streak. Parameters as in Fig. 15. 
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FIG. 24: (Color online) Growth rate, A, of the structurally unstable streaks in a Blasius boundary layer as a function of 
perturbation structure obliqueness, O = tan _1 (m/fc). The displacement thickness is Si — 1.72 and m is the wavenumber of the 
streak. Maximum growth rate occurs for m = 3 which corresponds to streak spacing Az = 2.4<5i or 50y + wall units consistent 
with observations. Also shown are the growth rates for m = 1,2,4. The channel width is L z = 4-7T, the Reynolds number is 
R — 400, the STM excitation parameter is / = 10 and the channel height is L y — 7. 



